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Magnetic Field Control of Exchange and Noise Immunity in Double Quantum Dots 
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We employ density functional calculated eigenstates as a basis for exact diagonalization studies of 
semiconductor double quantum dots, with two electrons, through the transition from the symmetric 
bias regime to the regime where both electrons occupy the same dot. We calculate the singlet-triplet 
splitting J{e) as a function of bias detuning e and explain its functional shape with a simple, double 
anti-crossing model. A voltage noise suppression "sweet spot," where dj{e)/de = with nonzero 
J(e), is predicted and shown to be tunable with a magnetic field B. 
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The goals of computation and information processing 
at the quantum level have stimulated the efforts of many 
researchers to coherently manipulate a variety of elemen- 
tary quantum systems. The scope of these candidate 
systems is wide |1|J. Advanced fabrication technology 
and inherent scalability, however, make semiconductor 
systems especially promising. The investigation of these 
systems has recently produced some auspicious results 
[2, 05 S H Q- The goal of these particular studies is 
to coherently manipulate and probe the spin and charge 
state of a small (typically electron number N=l or 2) 
system using: (i) time-varying electric fields from pulsed 
gates; (ii) charge sensors from nearby quantum point con- 
tacts (QPCs) 7, 8]; and (iii) externally applied magnetic 
fields B. Attention has focused recently on the regime 
adjacent to the degeneracy line between the double dot 
charge states (TVl, Nr) = (1,1) and (0,2) 2]. Here Nl 
and Nji denote the electron numbers on the left and right 
dots. 

In Ref. ^ the lateral gates confining the double dot 
were pulsed to produce a controllable "detuning" e of 
the potential {e is the potential difference between left 
and right gates measured from the degeneracy point of 
(1,1) and (0,2)), in order to first prepare two electrons 
in a singlet state in one (say, the right) dot, separate 
them into the two dots, and then recombine them in the 
right dot, i.e., (0,2) -^ (1, 1) -^ (0,2). For the employed 
gate voltages and dot level spacings, the recombination 
was suppressed by Pauli blocking in the case where the 
(1,1) electron is in a triplet 9]. The singlet-triplet split- 
ting, or exchange coupling J{s)^ the spin phase coherence 
time and the damping of Rabi oscillations between singlet 
and triplet were all thereby measured as functions of e at 
the separation point, the inter-dot tunnel coupling t, and 
magnetic field B. As £ ^ 0~, J{£) exhibited a rapid rise 
such that dJ/de increased with J. The effect of large 
dJ/de is to enhance sensitivity to voltage noise. The 
damping of Rabi oscillations, whose frequency is deter- 
mined by J(6:), appeared to increase with frequency, pre- 
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FIG. 1: (a) Self-consistent effective 2D potential profile at 
2DEG level, gate pattern superimposed; (b) Lowest two sin- 
glet and triplet CI energies and difference of lowest triplet and 
singlet J (triangles, right axis) versus detuning s, upper panel 
of (b) shows density of four states at £ = — 18my(dots indi- 
cated) with 5*0 and To delocalized and 5*1 and Ti mostly in 
right dot; (c) J{e) for two (singlet) tunnel coupling strengths 
t (triangles and boxes) and lowest singlet and triplet occu- 
pancies of left dot (right axis) for each t. 



sumably due to dephasing caused by this voltage noise. 

In this Letter we identify ranges of parameters where a 
noise-immunity "sweet spot" in exchange can be found, 
that is, where exchange is present (|J(£)| > 0) but 
the system is insensitive to first order electrical noise 
{dJ/de = 0). The sweet spot is identified in configuration 
interaction (CI) calculations M^ for the electronic struc- 
ture of the double quantum dot with A^ = 2 in Ref. [3 • 
The CI calculation employs basis states computed with 
density functional theory (DFT) pjj] to obtain full geo- 
metric fidelity to the experimental structure. We calcu- 
late J{s) from the symmetric limit at the center of the 
(1,1) honeycomb in the stability diagram well into the 
(0,2) regime with both electrons on a single dot. Analysis 
within a Hartree-Fock (HF), double anti-crossing model 
allows us to deduce simple expressions in the control pa- 
rameters by which the noise-immune regime can be ac- 
cessed I2J. 

DFT calculations for lateral heterostructures have 
been described extensively in the literature [ll| . We 
correct the DFT single particle energies of the N=2 dou- 



ble dot to avoid double counting of the Coulomb inter- 
action which is diagonalized in the basis of Kohn-Sham 
states 0i 0] . The advantage of this method is that the 
basis itself varies with gate voltages and B and thereby 
captures much of the evolving structure, much as a nat- 
ural basis does in quantum chemistry [13, [l3 • The basis 
states are states of the full double dot and so no artifi- 
cial tunneling coefficient needs to be incorporated |l5J . 
Finally, the Coulomb matrix elements automatically in- 
clude screening by the electrical environment and can 
be calculated very efficiently with the kernel of Poisson's 
equation, which is a natural byproduct of the DFT cal- 
culation. 

The self-consistent potential profile from the DFT cal- 
culation, with the device gate pattern superimposed, is 
shown in Fig. 1(a). Here, the bias is approximately 
symmetric and the potential minima of the two dots 
(~ 5 meV below the Fermi surface) are nearly equal. Fig- 
ure 1(b) displays the Cl-calculated lowest two singlet (S) 
and triplet (T) eigenstates as a function of e. The singlet 
and triplet ground states each anti-cross with their cor- 
responding first excited states (the calculation preserves 
total spin). Also plotted is J{e) = E^ — E^^ where E^ 
and E^ are the triplet and singlet ground states, respec- 
tively. The nature of the four anti-crossing states for 
e ^ —ISmV is shown by the total (2D) density plots on 
the figure. Here, near the anti-crossing, the two excited 
states on the (1,1) side each have their densities concen- 
trated in one dot, that is, they are the states that become 
S and T (0,2) ground states for £ > |1S]. Fig. 1(c) ex- 
hibits J{e) for two inter-dot tunnel-coupling strengths. 
Also shown are the S (red) and T (blue) total occupancies 
of the left dot versus e. The calculated results exhibit a 
coupling-dependent rapid increase of J(6:), as experimen- 
tally observed. Additionally, in the (0,2) region (e > 0), 
where the e dependence of J was not experimentally ex- 
plored, the calculation predicts a saturation. 

The structure of the numerical results in Fig. 1(b) sug- 
gest that we can explain the shape of J{£) in terms of two 
overlapping anti-crossings. These can be determined in 
their ideal form within a simple Hartree-Fock (HF) ap- 
proximation where the four "bare" states are written [I3 
as ^0 = {LI ^ Rl i), Si = {Rl T Rl i) (singlets) and 
To = (LI T i^l T), Ti = {Rl T R2 T) (triplets). The 
single particle energies are written e^^i, e^i and eR2 (we 
don't use 6^2)- We include the electrostatic potentials of 
the dot potential minima, eC(j)L and eC(j)R in the bare 
level energies (here C is a gate-to-dot lever arm, assumed 
constant) hence e^i = e/^i + eC^L, where e^j^^ denotes the 
energy measured from the dot bottom. Similar expres- 
sions hold for the other level energies. Within a HF type 
description, the bare eigenfunctions are retained and the 
energies are shifted by direct and exchange Coulomb ma- 
trix elements. Thus we write: E^ = e^i + €ri + Vinter^ 

2ei?l + Vintra, E^ = 6^1 + Cri + Vinter - Vf^^^^ and 

eRi^eR2^ Vintra " Vf^^^^ . For simplicity we here ig- 
nore the state-dependence of the inter- and intra- dot ma- 
trix elements. We also include the two (inter- and intra- 
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FIG. 2: (a) Schematic of double anti-crossing near (1,1) 
to (0,2) transition. Singlets (solid) anti-cross with tn, 
triplet (dashed) with ti2. Singlet-triplet splitting of ground 
(A5t(1, 1)) and excited {Ast{0,2)) states in (1,1) region in- 
dicated. Top inset: schematic of the configurations of the 
four anti-crossing states, (b) Fit of calculated J(s) to anti- 
crossing equation 121 (dashed line) for t = 0.17 meV case in 
Fig. 1. Triangles reproduce points in Fig. 1(c). 



dot) exchange matrix elements only in the triplet terms. 
The anti-crossing of the singlets results from the tunnel 
coupling of the two single particle ground states which 
we denote tu However, crucially, the triplets anti-cross 
with tunnel coupling ti2, since the Rl state is blocked 
by the Pauli principle. Due to the humped shape of the 
barrier, tunneling from LI to the higher R2 is stronger 
{ti2 > tu) and this is evident in Fig. 1(b). The two 
branches of the singlet and triplet are found by diagonal- 
izing the standard 2x2 determinants: 



Es± =Es ± 0.b^y{Ce)^^4t 
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Et± =Et ± 0.5^(5 -C£)2 + 4t 



(1) 
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where Et = {E^ ^ E^)/2 and Es = (£^g ^ E^)/^- As 
noted above, e is by definition zero at the singlet anti- 
crossing. In terms of the HF energies: Ce = eC{(j)L — 

(t>R) - ^Rl - CLl - Vinter + Vintra = 0. AlsO, S is the 

detuning (in units of energy, i.e. multiplied by the lever 
arm C) at which the triplet anti-crosses, relative to that 
of the singlet. It is given by S = e?^2 - 4i - "^in^'ter + ^i^ftra • 
The exchange splitting is the difference between the two 



lower branches inQ] 



J(e) = ET-Es-Q.^\{Ce - Sf + At\^+O.^J {Ce) 



(2) 
note that Et — Es is independent of e. 

Figure 2(a) shows the energy for the few levels, 
Es±-,Et±' Away and to the left of the anti-crossings, 
the gap between the ground states S and T is E^ — Eg = 
A5't(1, 1) = —yfnter H • The gap between the (bare) ex- 
cited states is E^-El = Ast{0, 2) = eR2 - em - Vf^^^^, 
i.e. it depends on a single particle level spacing in dot 
R. Clearly, from Fig. 2(b), as J increases near £ = so 
too does dJ/de. For spin manipulation with noise immu- 
nity it would be advantageous to find a regime where J 
was appreciable but dJ j de was not. This turns out to be 
possible. In Fig. 2(b) we also fit J{e) from Fig. 1(c) (i.e. 
the full CI results, specifically the curve with the smaller 
dot-dot coupling) with Eq. [21 The fit is good as long as 
points far into the (0,2) region are excluded jl9j . 

In what follows we show that the expression Eq. [2 
has exactly one extremum, dJ/de = 0, except where 
^11 = ^12- We then show in what parameter region the 
extremum is a minimum and how the parameters can be 
modulated by magnetic field B to accentuate that mini- 
mum. 

First, to show that J^e) has a single minimum, we take 
the derivative of Eq. [21 
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and observe that the two terms are sigmoidal curves 
which (for tn ^ ^12) must intersect in one point, specif- 
ically Em = ^/(l ~ (^12/^11))- A second derivative test 
shows that for ti2 > tn (the usual case), the extremum 
is a minimum. Note also that, assuming ti2 > tn, 
6 > ^ Em < ^ (the minimum is in the (1,1) zone) 
and ^ < =^ ^m > 0. 

Evaluating the depth D of the minimum of J{£) 
yields D = J{em) - J{e -^ -00) = 0.5^(1 - 
Vl + 4|ti2-tii|V|(5P ^ |ti2 - hi\y\S\. The fun CI 
calculations, which include the realistic inter-dot bar- 
rier, allow us to estimate the tunnel coefficients (com- 
pare Fig. 1(b) where the singlet and triplet anti-crossings 
have gaps of 2tii and 2ti2, respectively). Here, and in 
comparable parameter ranges, the tunnel coefficient ratio 
^12/^11 ^ 1-3. From its definition, 6 is the level spacing 
of dot R corrected by inter- and intra-dot exchange ener- 
gies and, at B=0, is of order 1.2 meV from calculations 
for the device in Ref. 2]. Note that S forms the major 
portion of A5't(0,2) (cf. Fig. 2(a)) and that reducing 
A5't(0, 2) while keeping Ast{^^ 1) relatively fixed causes 
the two anti-crossings to move closer to one another. For 
non-zero 5 in a circular parabolic potential, the lowest 
branch of the 2 state converges toward the 1 state and 
the exchange term (Hund's coupling) can induce a tran- 
sition in the N = 2 single dot to a triplet ground state; 
Fig. 4(a) and Ref. 20]. Simultaneously, B induces a 
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FIG. 3: dJ/de (solid, left scale) and J{e), from Eqs. 3 and 
2, resp., for tunnel and S parameters in Fig. 2(b). Note, 
offset energies, Es Et taken as zero throughout figure, so 
offsets of J arbitrary. Expanded view in (b) shows minimum 
{dJ{em)/de = 0) and depth of minimum. As £ ^ —00, J 
saturates above the figure. For same tunnel coefficients, dJ/de 
(c) and J{e) (d), for a range of values of 5. Minimum of J{e) 
(i.e. Sm) switches sign to positive when 6 changes sign to 
negative. 



singlet to triplet transition in the (1,1) ground state by 
enhancing the wavefunction overlap in the saddle point 
J^'>21,'>22i^ but the inffuence on the (0,2) single particle 
levels is greater and the result is a drastic decrease in 5, 
which can even become negative. 

Figure 4(b) shows J{e) computed (full CI) for a variety 
of B values and for t = 0.46 meV. The minima increas- 
ingly deepen as a function of B and move to larger e. 
Note that for this large value of tunnel coupling strength 
the minimum for small B (e.g. the B = 1.3 T curve) is 
not observed at all. In this case the tunnel coupling is 
comparable in scale to the level spacing and when 5 is 
large {B is small) mixing of higher levels apparently in- 
validates the simple, two- level HE picture which we have 
used. The point is, however, that t and B can be varied 
to change the order of singlet and triplet anti-crossings 
{5 goes from > to < 0). This shifts the minimum of 
J{e) from (1,1) to (0,2) and, for (5 = the minimum can 
become very deep. 

In summary, we have demonstrated that device param- 
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FIG. 4: Magnetic field dependence of J, from full CI calcu- 
lations, for center of (1,1) honeycomb cell (i.e. balanced gate 
voltages) triangles, and for strong (0,2) regime where both 
electrons occupy right dot, boxes. Note carefully difference 
between left and right scales. Single dot requires greater B to 
drive S-T transition, but \dJ/dB\ is much greater in (0,2). (b) 
Full CI calculation for J(s) for various B for strong coupling 
case in Fig. 1(c). Effect of B is principally to reduce single 
particle level spacing and thereby reduce Ast{0, 2) (compare 
Fig. 2(a)). This moves triplet anti-crossing to lower e, i.e. 
reduces 6, simultaneously deepening the minimum of J{e). 
When triplet anti-crosses before the singlet (^ < 0, B = 1.55 T 
and greater), the minimum of J{s) occurs in the (0,2) zone 
(s>0). 



eters can be tuned near the (1,1) to (0,2) crossover to 
generate a relatively noise- immune sweet spot where the 
exchange J{e) exhibits a minimum. The depth of the 
minimum, ~ \ti2 — tiip/|^|, and its detuning value, 
^m = ^/(l — (^12/^11)), can be modulated with magnetic 
field and tunnel barrier height. The existence of this 
noise-immune point for spin interactions is a promising 
development for quantum dot quantum computation. 
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